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Abstract. The NASA Discovery mission EPOXI, utilizing the Deep Impact flyby spacecraft, 
comprises two phases: EPOCh (Extrasolar Planet Observation and Characterization) and DIXI 
(Deep Impact eXtended Investigation). With EPOCh, we use the 30-cm high resolution visible 
imager to obtain ultraprecise photometric light curves of known transiting planet systems. We 
■ will analyze these data for evidence of additional planets, via transit timing variations or transits; 

^~~^ for planetary moons or rings; for detection of secondary eclipses and the constraint of geometric 

^ planetary albedos; and for refinement of the system parameters. Over a period of four months, 

^^ I EPOCh observed four known transiting planet systems, with each system observed continuously 

^12 ' for several weeks. Here we present an overview of EPOCh, including the spacecraft and science 

—^l , goals, and preliminary photometry results. 

r — , Keywords, planetary systems, space vehicles, methods: data analysis, techniques: image pro- 
{— ^ ■ cessing 
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• I— I . 1. Introduction 

rS The EPOXI mission, led by PI Michael A'Hearn and Deputy PI Drake Doming, is a 

d NASA Discovery Program mission of opportunity, designed to be a low cost, low risk, 

scientifically focused investigation. The mission comprises of two distinct science projects: 
Extrasolar Planet Observation and Characterization (EPOCh) and the Deep Impact 
extended Investigation (DIXI) . Both investigations make complementary use of the Deep 
Impact flyby spacecraft, an existing and flight-proven space asset that was previously used 
to observe comet Tempel 1. The EPOCh project is currently observing known transiting 
planet systems to obtain ultraprecise photometric light curves and is the focus of this 
contribution. 

The aim of the EPOCh project is to observe a carefully selected set of known, bright, 
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Figure 1. Left: The defocused point spread function of the high-resolution imager. Right: The 
pointing jitter for three days of HAT-P-4 data. The position of the center of the PSF in the 
128x128 sub-array is plotted, where points within 10 pixels of the edge of the chip have been 
discarded. The color scale indicates the passage of the star over the CCD with time over the 72 
hour period. 

transiting planet systems for several weeks at a time, obtaining light curves with very 
high precision, phase coverage and cadence. Using these data we can investigate a variety 
of science goals, including: refining system parameters; searching for additional bodies in 
the system via transit timing variations in the transits of the known planet (Agol et al. 
2005; Holman & Murray 2005) or detections of additional transits (CroU et al. 2007a, b); 
searching for reflected light at secondary eclipse and constraining the planetary albedo 
(Rowe et al. and Matthews et al., these proceedings); and searching for the signatures 
of moons (Doyle & Deeg, 2004, Brown et al. 2001) or rings (Barnes & Fortney, 2004) 
associated with the known planet. Another aspect of the EPOCh project was the char- 
acterization of Earth as an extrasolar planet (Deming et al. 2007); only the transiting 
planet observations and results arc discussed here. 



2. Observations 

We obtained the observations for EPOCh using the high-resolution imager (HRI) on- 
board the spacecraft. A telescope with a 30-cm aperture illuminates a clear filter with 
a bandpass of 350-950 nm. The light is imaged in 50 second integrations by a Ikxlk 
CCD with a pixel scale of 0.4" pixel" ^. Due to the limited amount of on-board memory 
(230MB), and our desire to maximize the time coverage between data downlinks from the 
spacecraft, only a 128x 128 sub-array of the CCD was initially used for the observations, 
resulting in a field of view of 51"x51". The left panel of Figure [1] shows the highly 
defocused point spread function (PSF); the advantages for high precision photometry 
include increased exposure time before saturation and the reduced effect of inter- and 
intra-pixel sensitivity variations. 

The fixed integration time of 50 seconds effectively limited the magnitude range in 
which our targets could lie to 9 > y > 13; our final target list and observing sched- 
ule are outlined in Table [T] These targets were selected to permit a range of scientific 
investigations. 

HAT-P-4 (Kovacs et al. 2007) is an interesting system due to its low planetary density. 
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Figure 2. The full time series of HAT-P-4 photometry, spanning 21 days. The gaps in the time 
series are due to periods when data are being downloaded from the spacecraft and also when 
pointing errors resulted in the target moving out of the field of view. 

which is unexpected given the high stellar metallicity and the presumed likelihood of the 
planet having a significant core. We observed seven transits of HAT-P-4, with six having 
better than 80% coverage for a duration three times the transit length centered on the 
time of transit, which is the baseline requirement. The full time series of HAT-P-4 is 
shown in Figure [2l 

XO-3 (Johns-KruU et al. 2008) is one of the most massive transiting planets known to 
date and lies in the relatively unpopulated region between high-mass planets and brown 
dwarfs. No transits of XO-3 were observed due to the spacecraft traversing perihelion 
during this time. This increased the on-board temperatures beyond the nominal operating 
range and the spacecraft entered safe mode. 

TrES-3 (O 'Donovan et al. 2007) has an extremely short period and is an excellent 
target for reflected light detection, however only 5 transits and 6 secondary eclipses of 
TrES-3 were able to be observed due to the time lost before the spacecraft was recovered 
from safe mode. 

The spacecraft then experienced a significant decrease in the telecommunications power, 
and for our next target, XO-2 (Burke et al. 2007), which is one of the few known transiting 
planets to be found in a visual binary star system, we were only able to download partial 
data for three transits. At this stage observations were suspended until the loss in power 
could be explained and observations on our penultimate target, TrES-2 (O'Donovan et 
al. 2006), were not obtained. TrES-2 is in the Kepler field of view (Borucki et al. 2003), 
and will be observed continuously from space for at least three years by that mission, 
starting in 2009. Contingent EPOCh observations of TrES-2 (see below) will provide an 
increased baseline of precise timing measurements with which to probe transit timing 
variations. 

The system was re-tested when the Earth-spacecraft distance had decreased from the 
maximum separation, and the power was found to have returned to full downlink capacity. 
We then resumed observations on our final target in the original schedule, GJ436 (Gillon 
et al. 2007) . This is an exciting prospect for further study due to the non-zero eccentricity 
of the planet; one explanation is the presence of the third body in the system, pumping 
up the eccentricity. Seven transits of GJ436 were observed to complete the initial schedule 
of observations. 

Besides the losses due to the spacecraft downtime, there were also significant losses due 
to pointing errors. Stable pointing is highly desirable for very high precision photometry, 
and we found that the pointing was insufficiently accurate, in that not only did the 
initial pointing at the target have a significant offset from the center of the CCD, but the 
subsequent pointing drifts resulted in the target either leaving the field of view entirely or 
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Table 1. EPOCh target list and observing schedule. 



Target 


V Magnitude 


Period (d) | Scheduled 


Transits observed 


HAT-P-4 


11.0 


3.06 


1/22/08-2/12/08 


6 


XO-3 


9.9 


3.19 


2/13/08-2/19/08 





TrES-3 


12.4 


1.31 


2/20/08-3/18/08^ 


5 


XO-2 


11.3 


2.62 


3/20/08-4/07/08 


0^ 


TrES-2 


11.4 


2.47 


4/08/08-4/29/08 





GJ436 


10.7 


2.64 


5/01/08-5/28/08^ 


7 


Contingent observations 






HAT-P-4 


11.0 


3.06 


6/29/08-7/07/08^ 


- 


TrES-2 


11.4 


2.47 


7/08/08-7/29/08-^ 


- 


WASP-3 


10.6 


1.84 


7/30/08-8/15/08 


- 


HAT-P-7 


10.9 


2.20 


8/16/08-8/31/08 


- 



Notes: 

The baseline requirement for a successful transit observation was 80% coverage. Three partial transits of 
XO-2 were observed. 
^ Earth observations were obtained on 3/19/08, 5/29/08 and 4/06/08. 

These targets are being returned to in order to fulfill mission baseline criteria. 

lingering near the edge, where it is difficuh to recover accurate photometry. An example 
of the pointing wander is shown in the right panel of Figure (TJ where the position of the 
center of the PSF on the CCD is plotted for over three days of observations (Rieber & 
Sharrow, in prep). Additionally, due to the delay in downloading the data, calculating 
the required pointing correction and uploading the new coordinates, several days to a 
week of observations would have been completed in the meantime without correction. 

Two strategies were employed to alleviate the pointing problem. The first was the 
inclusion of a 'pre-look', whereby the spacecraft is pointed at a new target, sufficient 
data obtained to deduce the pointing offset, and the spacecraft then returned to observe 
the old target while the pointing correction on the new target was calculated. In this way, 
no observing time was used inefficiently waiting for the upload of the pointing correction. 
The second alteration was the increase in the dimensions of the sub-array from 128x 128 
to 256x256 during the predicted transit and eclipse windows, to ensure photometry could 
be obtained at these critical times. This was put into place for the observations of GJ436 
and will be used in the contingent observations discussed below. 

The investigation was initially allocated slightly more than four months of observing 
time, two of which were effectively lost during safe mode and the following telecommu- 
nications power loss. In order to recoup the science losses incurred during this period 
reserve funds were made available to obtain an additional two months of observations 
beyond the original schedule. These contingent observations are included in Table[T] Our 
plan is to return to HAT-P-4 brieffy to complete the mission baseline of five secondary 
eclipses, and to observe TrES-2, one of the scheduled targets missed in the spacecraft 
downtime. The final two targets that we will observe are WASP-3 (PoUacco et al. 2008) 
and HAT-P-7 (Pal et al. 2008), both very good targets for reflected light detection. The 
latter is also in the Kepler field of view. 



3. Data reduction 

The raw data are converted to calibrated images in a standard fashion (bias- and dark- 
subtracted and flat-fielded) using the pre-existing Deep Impact data reduction pipeline 
(Klaasen et al. 2005). The extraction of the photometry is summarized here and is de- 
scribed in more detail in Ballard et al. (these proceedings). We use a drizzled PSF, 
oversampled by a factor of 10, to locate the stellar position; bilinear interpolation then 
results in positions estimated to within one hundredth of a pixel. PSF-fitting to produce 
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the final photometry is complicated by the underlying structure of the CCD, where the 
central two rows appear to be physically smaller, altering the shape of the PSF when 
it straddles these rows. We are currently analyzing the data using aperture photometry 
and scaling each of the rows by a constant factor that is individually determined for each 
target so as to minimize the rms of the time series. 

Initially we found the light curves contained a significant component of red noise, 
correlated with the star's position on the CCD. The flat fields used in the calibration 
pipeline were obtained in the lab when the spacecraft was on the ground, and the flat 
field structure of the CCD may have changed significantly during the flight. An on-board 
'stimulator' lamp comprised of a green LED mounted near the CCD was used in an 
attempt to create an updated flat fleld, however we have found much better results using 
a decorrelation of the data with position, based on a local 2D spline fit. This has largely 
removed the red noise component from the light curve; the data shown in Figure [2] are 
the decorrelated data. The next dominant noise source is the presence of a population 
of low-energy radiation events, which can be clearly seen in Figure [5] as a preponderance 
of outliers above the median. The noise in the decorrelated data now decreases with bin 
size as V^, where N is the number of points per bin, as expected for a Gaussian noise 
profile. 



4. Preliminary analysis and conclusions 

Our preliminary analysis has thus far focused on the HAT-P-4 dataset; the initial 
investigations and further work are briefly described below. 

4.1. Refinement of system parameters 

We have used the analytic equations describing transiting light curves of Mandel & Agol 
(2002) and the Monte Carlo Markov chain (MCMC) method outlined for radial velocity 
measurements by Ford (2005) and adapted to transit light curve measurements by, e.g., 
Holman et al. (2006) to analyze the light curve of HAT-P-4. We used the four-parameter 
non- linear limb-darkening law of Claret (2004), with the limb-darkening coeHicients de- 
rived from a weighted fit across the instrument bandpass, interpolated across the closest 
theoretical stellar modelij] (Kurucz 1994, Kurucz 2005). The phase-folded binned light 
curve is shown in Figure [3] with the best-fit solution over-plotted. Further investigation 
of the correlation of the errors in the parameters is required, however our results are thus 
far consistent with those presented in the discovery paper (Kovacs et al. 2007). 

4.2. Reflected light analysis 

We have analyzed the data around secondary eclipse in order to search for the signature 
of reflected light and to put constraints on the geometric planetary albedo. The MOST 
satellite has been used very successfully to put tight constraints on this parameter for 
HD 209458b and HD 189733b (Rowe et al. 2006, see also the Rowe et al. contribution to 
these proceedings). The upper panel of Figure [3] shows the phased HAT-P-4 data around 
the predicted times of secondary eclipse (set to a phase of 0.0), binned by a factor of 80. 
The best-fit eclipse model gives an upper limit on the geometric albedo of 0.346, below 
that of Jupiter, 0.6. Further work allowing the predicted time of eclipse to vary owing to 
the range in eccentricity allowed by the radial velocity measurements is still required. 



t http://kurucz.cfa.harvard.edu/| 
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Figure 3. Phase-folded HAT-P-4 data, overplotted with the best-fit analytic solution from Man- 
del & Agol (2002), with fixed limb-darkening coefficients derived for the instrument bandpass. 
The data are binned by a factor of 20 and the error bars are calculated using the out-of-tran- 
sit standard deviation. The residuals to the solution are shown underneath, with a standard 
deviation of 0.4 mmag. 

4.3. Transit timing variations 

The timestamps for the HAT-P-4 data have been corrected to the barycentric Julian 
date (BJD), however the relative offset between the clock on-board the spacecraft and 
UT on Earth is yet to be applied. As a result, we can conduct an independent period 
search and transit timing analysis within the EPOCh dataset, however at this stage we 
cannot combine these results with those of ground-based observations. We find a period of 
3. 056370±0. 000094 days, and Figure [5] shows the residuals of the individual transit times 
to this period. When the timing offset has been applied, we will analyze the long-term 
period and transit timing behaviors. 



References 

Agol, E., Steffen, J., Sari, R., & Clarkson, W. 2005, MNRAS, 359, 567 

Barnes, J. W., & Fortney, J. J. 2005, ApJ, 616, 1193 

Borucki, W. J., et al. 2003, ProcSPIE, 4854, 129 

Brown, T. M., Charbonneau, D., Gilliland, R. L., Noyes, R. W., & Burrows, A. 2001, ApJ, 552, 

699 
Burke, C. J., et al. 2007, ApJ, 671, 2115 
Claret, A. 2000, A&A, 363, 1081 
CroU, B., et al. 2007a, ApJ, 658, 1328 
CroU, B., et al. 2007b, ApJ, 671, 2129 

Deming, D., et al. 2007, AAS/Division for Planetary Sciences Meeting Abstracts, 39, 2202 
Doyle, L. R., & Deeg, H. J. 2004, lAUS, 213, 80 
Ford, E. 2005, AJ, 129, 1706 
Gillon, M., et al. 2007, A&A, 472, L13 
Holman, M. J., & Murray, N. W. 2005, Science, 307, 1288 
Holman, M. J., et al. 2006, ApJ, 652, 1715 
Johns-KruU, C. M., et al. 2008, ApJ, 677, 657 



The NASA EPOXI Mission 



125 



1.0010 
g 1.0005 
^ 1.0000 



0.9990 
-0 



IftMWii+V/ 



.10 



-0.05 



0.00 
Phase 



HAT-P-4 Reflected Light 
J 1 1 I , 1 1 

0.05 



0.10 





ZU 


:' ' \ ' ' ' ' ' ' 


, , , 1 


: 


T3 

rr 


15 


- \ 


/ 0.00007128 


— 


in 




\ 


/ 


_ 


Jl 




\ 


/ 


- 


O 


10 


— ^^ 


/ 


— 


n 


\ / 




_ 












Q 


5 


r \. y 




- 







- 1 ^^^^-^ ^.^^1 


1 


z 



-0.0002 



0.0000 
Eclipse depth 



0.0002 



Figure 4. Searching for tlie secondary eclipse of HAT-P-4. Upper. Tlie data around secondary 
eclipse (set to a phase of 0.0) are shown here, binned by a factor of 80. The solid line shows 
a model with a secondary eclipse depth inconsistent with the data at a 3(t significance, giving 
an upper limit on the planet-to-star flux ratio of 7.1 x 10~^. Lower. The change in chi-squared 
from the best-fit value, with the vertical dashed line delineating the allowed upper limit on the 
depth. 
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Figure 5. Residuals of the measured center of each individual transit of HAT-P-4 to the 

best-fit period. 
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